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A METHOD OF APPROXIMATION. 



By S. A. COREY, Hiteman. Iowa. 



The following method of approximation may or may not be new, but as I 
believe it to be of practical importance I wish to call attention to it, and to point 
out that to secure the same degree of accuracy, it involves, at least in some cases, 
less labor than does the method of mechanical quadrature. It can also be used 
in certain cases where the method of mechanical quadrature fails, because the 
finite differences used in the latter method cannot always be found between the 
proper limits. On account of its rapid convergence it can also be used when 
other common methods fail. 

We have the formula,* 

A.)=A0)+^{/W+/(0)+2[/(^) + /'(|i) + ...4/(^)]} 

+(-D n ^9o! C/2 " (X)_/2 " ( ° )] + (1) ' 

(2?!, Ba, i? 3 , etc., being Bernoulli's numbers, \, -fa, T V> to» A> etc 
For brevity, (1) may be written, 



ftz)=F(z,m)+-^+-^+-=*:+ etc. ...(2), 

in which s 2 , s 4 , s e , etc., are independent of m. It is evident that m may be con- 
sidered a convergence factor which makes F(x, m) approach /(a:) as m approaches 
infinity, F and/ becoming equal as m becomes infinite. But as m increases, the 
labor of evaluating F(x, m) also increases, so that, in practice, it is found desir- 
able to obtain one or more of the higher derivatives involved in s s , s i , s e , etc., 
in order to permit the use of a smaller value of m, and yet to attain the required 
degree of accuracy. In certain cases it becomes exceedingly difficult to obtain 
the form and value of the higher derivatives involved in s s , s 4 , s 6 , etc., and 
chiefly for this reason, it becomes desirable to eliminate from (2) certain of these 
quantities, (s 8 , s 4 , s e , etc.), after a certain degree of approximation has been 
reached. To accomplish this elimination take more than one value of m and set 
down a separate equation for each value chosen, thus, 

*See Annals of Mathematics , Second Series, Vol. 5, No. 4, July, 1904. 
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f(x)=--F(x, m 3 )+-%+-^4+-^ + (5), 
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etc. (In selecting the values of m it is usually convenient to take the highest 
value such that one or more of the lesser values will be stib-nmltiples thereof.) 
It is evident that from (»"+l) such linear equations there may be elimi- 
nated r of the quantities s, say s 2i , Sj^+u, ..., s a( i +r _ 2 ) . There results, 
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...(6), 



where M 1} M g , ..., Jf( r +i) , represent the approximate value of/(«) obtained by 
(2) for m 1 , m 2 , ..., m (r+ 2) , respectively, by the use of all the terms preceding 

If all derivatives higher than the (2i-j-2r— l)th are zero, or if all such 
even-numbered derivatives are equal for x=x and x—0, (6) gives the exact value 
of f(x). In other cases the degree of accuracy attained will be no less than 
would have been obtained by taking the smallest value of m used in (6) and de- 
veloping f(x) by (2) so far as to include all the r terms eliminated in (6). 

As a special case of (6), let r=2, «=2. Here M=F(x, m)-\ — ^, and we 

get after reducing, 

m 1 6 (w 2 2 -w 3 8 )Jf 1 +m a 6 (m 3 8 — m*)M t +mf(mf — mf)M 9 _ 
•^ m, 6 (m 8 2 — m i 2 ) + m e 6 {m^— m?)-\-mg(mf — m 2 a ) 



*Determinants of this form may be readily evaluated. See Weld's Theory of 
Determinants, Articles 23 and 27. 
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It is clear that the more nearly all the STs approach f(x) , the more nearly 
does (6) give the exact value of /(#), and that this increase in approximation is 
obtained by sufficiently increasing the integers, to, i, and r, or any of them. 

The following points may be noted in passing : 

First. After choosing certain values of m, r, and i, f(x) is developed by 
(6) in linear terms of the (r+1) quantities (M l ,M. i , ..., M^ +1) ), each M being 
an approximate value of f(x). 

Second. If M—F(x, to) and if the values (but not the form) off(x) be 
given at each of the (m + 1) equidistant points (0, x/m, 2x/m, ..., x) for each of 
the (r+1) values of to chosen, the value of f(x) may be very closely approximat- 
ed, although in such case the law of its development would not be expressly 
stated, but it must be assumed that some such law exists in order that/(«) may 
be developable in terms of x. This makes the method of practical value in many 
scientific problems where it becomes necessary to find f(x) approximately, and 
when the only available data are a number of observations of the values of f(x) 
at intervals of x/m, the approximation being closer in this method than in that 
of mechanical quadrature, especially where the intervals (x/m) are few. 

A discussion of the cases where (6) does not give the value of f(x) is not 
here entered into, but the following considerations will be pertinent : 

First. As (6) is obtained from a number of equations, (3), (4), (5), 
etc., care should be taken that each of these equations actually does give a closer 
approximation of f(x), when the number of terms employed include those elim- 
inated in deriving (6), than does M. 

Second. As (6), (2), and (1) are all dependent on the sum of a number 
of series developed by Stirling's formula, it becomes necessary that each of these 
underlying Stirling's series should give a true and convergent development of / 
between any and all of the adjacent, equidistant points (x/m) at which the value 
of / is taken. 

The following example has been chosen because its development by form- 
ula (1) has been given in the Monthly,* Vol. XIII, No. 4, and because its solu- 
tion by mechanical quadrature has been given by Dr. G. W. Hill in the Analyst, 
Vol. II, p. 120, 1875.t 



Example. Evaluate ( 



4* xdx 



o sina;(l+.16cos 2 a:)5 

Taking r=2, «=2, to, =6, to 2 =3, to 3 =2, we get M 1 =1.657,626,355, 
M s =1.657,490,406, M 3 =1.657,013, 853. Substituting in (7), we get, 

*The value of the integral is incorrectly given in the Monthly as, 1.657,636,524, in- 
stead of 1.657,636,257. The error is due to the tact that 



r 2596*' -1 
™' L4200.12«.6!J 



is carried out as .000,000,009 instead of .000,000,276. 

tSee also Hill's Collected Works, Vol. I, p. 204. 



140 

fa - ff 0^.= f 4 " xdx — 1729M * - 2332 83T,+ 2332803^ 

A**) AV)-J q s i n:r (i + .it5cos 2 a;)3 ~ 211,680 

=1.657,636,33.. .(8). 

We may check the accuracy of the work of computation and determine the 

degree of approximation attained in (8) by finding the values of jT(-i-) and 

3r 
/' ( -Q-). This done we may take mj=6, m 3 — 4, m 2 =3. Substituting these 

values of m t , m 2 , rw 3 , in (7) we get an expression similar to (8) which gives a 
value of the definite integral coinciding with the value found in (8) for six deci- 
mal places, thus proving the accuracy of the work and showing the number of 
decimal places to which result found by (8) is correct. 

This result is not as accurate as that obtained in the April Monthly for 
the reason that the smallest value of m here used is 2. This is therefore a more 
accurate result than would have been obtained by taking m=2 and developing 
by (1) so far as to include the term involving B t , but a less accurate result than 
was obtained in the Monthly by taking m—6 and developing by (1) far enough 
to include the term involving B 3 . The above result is, however, as accurate as Dr. 
Hill's and was obtained with much less labor. To obtain a more accurate result 
than that given in the Monthly without finding any higher derivatives than are 
there given, take i=A, r=l, m 1 =:6, m 2 =3, and substitute in (6). By using nine 
decimals throughout the result is found to be 1.657,636,259. By using ten deci- 
mals throughout a result correct to nine or ten decimal places would have been 
obtained. 



APPROXIMATION OF THE GREATEST ROOT OF A CUBIC 
EQUATION WITH THREE REAL ROOTS. 



By CHARLES GILPIN, JR., Philadelphia, Pa. 



We are concerned with the "irreducible case," in which Cardan's formula 
is of no value for computation. By replacing x by —x if necessary, we need 
consider only the form 



r.3 — 



ax—b—0 (a and b positive) ...(1), 

the two lesser roots of which are negative and the greatest root, which equals 
the sum of the lesser ones, is positive. It can be shown* that the greatest root 
g lies between the limits -\/a and i/(4a/3). From equation (1) we obtain 

*g> J a by (2). Since 6 2 /4 — aS /27<0 in the irreducible case, x* — ax — b is posi- 
tive for #=J(4a/3), negative for x=4a. Editor. 



